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Abstract. We define and study in detail utraslow scaled Brownian motion (USBM) 
characterised by a time dependent diffusion coefficient of the form D(t) ~ 1/t. For 
unconfined motion the mean squared displacement (MSD) of USBM exhibits an 
ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled 
Brownian motion. In an harmonic potential the MSD of USBM does not saturate 
but asymptotically decays inverse-proportionally to time, reflecting the highly non¬ 
stationary character of the process. We show that the process is weakly non-ergodic in 
the sense that the time averaged MSD does not converge to the regular MSD even 
at long times, and for unconfined motion combines a linear lag time dependence 
with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms 
of the ergodicity breaking parameter. The USBM process is also shown to be ageing: 
observables of the system depend on the time gap between initiation of the test particle 
and start of the measurement of its motion. Our analytical results are shown to agree 
excellently with extensive computer simulations. 


1. Introduction 

In the wake of the development of modern particle tracking techniques strong deviations 
of the time dependence of the mean squared displacement (MSD) from the linear law 
(. x 2 (t )) ~ t derived by Einstein [TJ and Smoluchowski [2] have been observed in a variety 
of complex fluidic environments 0IH El El El- Typically, anomalous diffusion of the 
power-law form 

(x 2 (f))~f“ (1) 

is observed, where, depending on the value of the anomalous diffusion exponent a, 
we distinguish subdivision with 0 < a < 1 and superdiffusion with a > i mu- 
Accordingly, subdiffusion was observed in the cytoplasm of living cells nm mi , in 
artificially crowded liquids pa mu, and in structured or functionalised environments 
[QJ. Also superdiffusive motion was found in living cells [151 fl6]. 
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Recently, interest in ultraslow diffusion processes with the logarithmic form 

(x 2 (t)) ~ log 7 (f) ( 2 ) 

of the MSD with different values for the exponent 7 has been revived [ 6 ]. Ultraslow 
diffusion may be generated by periodically iterated maps m and observed for random 
walks on bundled structures [T 8 ] • A prototype model for ultraslow diffusion is provided 
by Sinai diffusion in quenched landscapes with random force held, for which 7 = 4 
[19, I2DI EH [22J. In the context of Sinai diffusion ultraslow continuous time random 
walks with super heavy-tailed waiting times with 7 > 0 [22], S3, E 3 J [ 2 S] were discussed. 
Ultraslow scaling of the MSD of the form (|2]) were obtained in aperiodic environments 
(variable 7) [2S] and vacancy induced motion (7 = 1) [2Tj . Moreover, it occurs in 
heterogeneous diffusion processes with exponentially varying diffusivity (7 = 2 ) [28]. 
or interacting many-body systems in low dimensional disordered environments with 
7 = 1/2 [29], the dynamics of the latter being governed by an ultraslow, ageing counting 
processes [30] . 

The logarithmic time dependence (J2]) with 7 = 1 of the MSD is also observed 
for the self diffusion of particles in free cooling granular gases with constant, sub¬ 
unity restitution coefficient in the homogeneous cooling state [HI]. Granular gases are 
rarefied granular systems, in which particles move along ballistic trajectories between 
instantaneous collisions [31]. They are common in Space, for instance, in protoplanetary 
discs, interstellar clouds and planetary rings [32]. At terrestrial conditions granular 
gases may be obtained by placing granular matter into containers with vibrating [S3] 
or rotating [31] walls. If no net external forces (gravitation, etc.) are acting on the 
granular system, the motion of granular particles gradually slows down due to dissipative 
collisions between them [31] . This microgravity condition can be achieved, inter alia, 
with parabolic airplane flights or satellites [35], 36] 07] or by the use of diamagnetic 
levitation |38j. We note that in very dense two-dimensional lattice gas systems, ultraslow 
diffusion emerges, as well [39]. 

Figure [I] shows the crossover from the ballistic to the ultraslow form ([2]) of the MSD 
of a granular gas with constant restitution coefficient e = 0.8 in the homogeneous cooling 
state. Haff’s law demonstrates that the kinetic temperature of such a free granular gas 
with constant restitution coefficient decays inverse-proportionally with time, T(t) — l/t 
[4D]. For the effective self diffusion of the gas particles—mediated by particle-particle 
collisions—this property translates into the time dependent diffusion coefficient D{t) — 
l/t [H, 021 33 ] • We note that a diffusivity of the form D(t ) = D 0 + Di/t with a 
component decaying inverse-proportionally with time was used in the modelling of the 
motion of molecules in porous environments [44] as well as of water diffusion in brain 
tissue measured by magnetic resonance imaging m- 

Here we study in detail the process of ultraslow scaled Brownian motion (USBM) 
with time dependent diffusion coefficient D(t) — l/t. Starting from the Langevin 
equation for USBM and a summary of the simulations procedure we present analytical 
and numerical results for the MSD and the time averaged MSD for the cases of 
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Figure 1 . Time dependence of the ensemble averaged MSD obtained from 

event driven molecular dynamics simulations of three-dimensional force-free granular 
gases [4H| . At short times the particles follow ballistic trajectories, while for longer 
times the ensemble averaged MSD has a logarithmic time dependence. The inset 
focuses on the logarithmic long time behaviour. 


unconfined (Section [2]) and confined (Section [3]) motion. We analyse in detail the 
disparity between the ensemble and time averaged MSD and quantify the statistical 
scatter of the amplitude of the time averaged MSD of individual realisations of the 
USBM process. Moreover we study the ageing properties of USBM, that is, the explicit 
dependence of the physical observables on the time difference between the initiation of 
the system and the start of the observation. In Section [4] we present our Conclusions. 
In the Appendix we present details of the calculation of higher order moments and the 
ergodicity breaking parameter. 

2. Unconfined ultraslow scaled Brownian motion 

2.1. Overdamped Langevin equation for ultraslow scaled Brownian motion 

Anomalous diffusion processes with power-law form ([Tj) of the MSD are often modelled in 
terms of scaled Brownian motion (SBM) characterised by an explicitly time dependent 
diffusivity of the power-law form D(t) — t a ~ l with 0 < a < 2, see, for instance, 
references PI SUSHI ESI EniEU as well as the study by Saxton (52] and further references 
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Figure 2. Schematic of the motion of a Brownian particle in a bath with decreasing 
temperature T(t) — t 2a ~ 2 for 0 < a < 1. The diffusion coefficient of the Brownian 
particle decays with time as D(t ) ~ f Q_1 . USBM corresponds to the case a = 0, while 
standard SBM is strictly limited to 0 < a < 2 [53] . 


therein. In SBM this form of D{t) is combined with the regular Langevin equation [53] 

^ = v/2 W) X CM, (3) 

in which ((t) represents white Gaussian noise with the normalised covariance 

(c(t 1 )at 2 )) = 5(t 1 -t 2 ) ( 4 ) 

and zero mean (</(f)) = 0. While for a system connected to a thermal reservoir a 
description in terms of a time dependent temperature underlying SBM is unphysical 
|54J, time dependent diffusion coefficients appear naturally in systems that are open or 
dissipate energy into other degrees of freedom such as the granular gases discussed above, 
see the schematic in figure [2} In fact, granular gases with a viscoelastic, relative particle 
speed-dependent restitution coefficient correspond to SBM with a = 1/6 [31] 143]. 
Diffusion in media with explicitly time dependent temperature can, for instance, also 
be observed in snow melt dynamics [551 [56] . 

A diffusion equation with a time dependent diffusivity proportional to t 2 was 
originally introduced by Batchelor EH to describe the anomalous Richardson relative 
diffusion [53] in turbulent atmospheric systems. SBM with diffusivity D[t) ~ t a ~ l was 
studied extensively during the last few years [59 ] 160 ] l6l ] 154 ] 621) . In particular, the weakly 
non-ergodic disparity between ensemble and time averages in SBM as well as its ageing 
behaviour were analysed [hU, SB EH, 62], see also below. Processes with both time and 
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position dependent diffusion coefficients were also reported [63j. SBM is a Markovian 
process with stationary increments £(t), however, it is rendered non-stationary by the 
time dependence of the coefficient D[t). SBM is therefore fundamentally different 
[61(51! from seemingly similar processes such as fractional Brownian motion or fractional 
Langevin equation motion (63]. 

Following the motivation from our studies of granular gases with constant 
restitution coefficient [U] we here consider USBM with the time dependent diffusion 
coefficient 

Do 

1 + t/r 0 ' 

The time scale To defines the characteristic time beyond which the long time scaling 
D(t) D 0 r 0 /t sets in. We here introduce r 0 to avoid a divergence of D{t ) at t = 0. The 
case ([5]) is explicitly excluded in the allowed range for the scaling exponent a in SBM and, 
as we will see, constitutes a new class of stochastic processes. In the following we solve 
the overdamped Langevin equation (|3]) with the time dependent diffusion coefficient (J5]) 
analytically and perform extensive computer simulations of the corresponding finite- 
difference analogue of the Langevin equation. In this procedure, at each time step the 
increment of the particle position takes on the value 

x i+ i - Xi = \j2D{i){W i+ i - Wi), i = 0 , 1 , 2 ,..., (6) 



where W i+1 — W t is the increment of the standard Wiener process and D{i) is the value 
of the time dependent diffusivity ([H]) at the time instant i. We simulated N = 10 3 
independent particles (runs) with the parameters To = 1 and D 0 = 1/2 in all graphs 
presented below. 


2.2. Ensemble and time averaged mean squared displacements 


From direct integration of the Langevin equation (|3]) with the time dependent diffusivity 
(J5| we find the ultraslow, logarithmic growth 

(x 2 (t)) =11 \/ D{t')D (i t") (C(t')C (t")) dt"dt' 

Jo Jo 

= 2D 0 t 0 log ^1 + (7) 


of the ensemble averaged MSD. USBM therefore reproduces the asymptotic behaviour 
of the MSD for granular gases in the homogeneous cooling state and with constant 
restitution coefficient [43], as shown in figure [I] 

In addition to the ensemble averaged MSD ( x 2 (t )) of the particle motion, it is often 
useful to compute the time averaged MSD 


r t -a r 


h 2 (A) = 


x{t' + A) — x(t') 


dt!. 


( 8 ) 


t ~ A Jo 

Here, the lag time A defines the width of the averaging window slid over the time series 
x(t) of the particle position of overall length t (the measurement time). Time averages 
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of the form ([8]) are often used in experiments and large scale simulations studies based 
on single particle tracking approaches, in which typically few but long trajectories are 
available uni hd Eg. The careful analysis of the time averaged MSD (|8]) provides 
additional important information on the studied process as compared to the ensemble 
averaged MSD ( x 2 (t )), see, for instance, the analyses in references [TT1155]. Often one 
takes the additional average over N individual particle traces 5 2 (A), 

(W(&)) =-'£5?{A). (9) 

i= 1 

For ergodic processes^ such as Brownian motion, fractional Brownian motion, and 
fractional Langevin equation motion the time averaged MSD converges to the ensemble 
averaged MSD in the limit of sufficiently long times, lim^ 0O 5 2 (A) = (x 2 (A)) [6]. This 
property is due to the stationarity of the increments of these processes |66j|. The ergodic 
behaviour lim^oo 5 2 (A) = (x 2 (A)) of these processes holds for unconfined motion when 
the system is in fact out-of-equilibrium, an advantage of the particular definition ([8]). 
Moreover, ergodic systems fulfil the equivalence 


(i 2 (A)) = h 2 (A)> (10) 

even at finite t [Oj. Systems in which we observe the disparity ^h 2 (A) j ■/- (x 2 (A)) and 

therefore also Hindoo <5 2 (A) ^ (x 2 (A)) are called weakly non-ergodic [T; '5j [6, Tj [67, 68] j§] 
To calculate the time averaged MSD @i for USBM we do not need to consider 
the mixed position autocorrelations in the definition of the time averaged MSD, as the 
expression in the angular brackets simplify as follows, 


5 2 (A) = 


t - A 
1 


rt~ A 


ct- A 


x{t' + A) — x{t') 


t - A 

This is due to the propert}|||] 


{x (if + A)) - (x (if)) 


dt' 


dtf. 


( 11 ) 


(x(t)x(t + A)) = (x 2 (t)). 


( 12 ) 


for stochastic processes whose increments are independent random variables. We thus 
find the exact form for the time averaged MSD of USBM, 



2DqTq 

t - A 


£(t) - £(A) - £(t - A) , 


(13) 


f We consider processes ergodic in the Boltzmann-Khinchin sense when the long time average of a 
physical observable converges to the associated time average. 

§ Note that also transiently non-ergodic behaviour may become relevant as it may mask intrinsic 
relaxation times when time averages are measured }12[ l69j . 

|| In contrast, this is not valid in the case of granular gases, where particles move ballistically in 
between instantaneous collisions [43j . or for processes driven by long-range correlated increments such 
as fractional Brownian motion or fractional Langevin equation motion [61 EH [70]. 
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where we introduced the auxiliary function 
£(t) = (t + r 0 ) log (l + — 

V To 


The time averaged MSD (13) thus crosses over from the limiting behaviour 


5 2 ( A) 


A 


2A)Toy log ( — 


(a 


(14) 


(15) 


at short lag times Tq <C A <C t combining a linear with a logarithmic A dependence, to 
the purely logarithmic law 

t +t 0 


(V(A)^) ~ 2D 0 t 0 log f - 


. . . (16) 
\t — A + To J 

at To C A « t We see that as the lag time A approaches the measurement time 
t, the time average MSD approaches the MSD (7), ^ 5 2 {t — > (x 2 (t)). The results 
of our simulations of the USBM process for both ensemble and time averaged MSDs 
agree very well with the above analytical results, as demonstrated in figure [3} In that 
plot the thin grey curves depict the simulations results for the time averaged MSD for 
individual trajectories. The amplitude spread between different trajectories is fairly 
small for A <C t and increases when the lag time A approaches the trace length t due 
to worsening statistics. 


2.3. Stochasticity of the time averaged mean squared displacement and ergodicity 
breaking parameter 

Even ergodic processes such as Brownian motion exhibit a certain degree of stochasticity 
of time averaged observables for shorter measurement times. The amplitude fluctuations 
at a given lag time A of the time averaged MSD as compared to the trajectory average 
0 is quantified in terms of the ergodicity breaking parameter m edi m m 


EB(A) = lim 

t—¥ OO 


<5 2 (A) 


- (5 2 (A)) 


y(A) 


= lim (£ 2 ) - 1, 

t—¥ OO 


(17) 


where in the second equality we introduced the relative deviation m 
£ = 


h 2 (A) 


5 2 ( A) 


(18) 


The necessary condition for ergodicity of a stochastic process is that the ergodicity 
breaking parameter vanishes in the limit of infinitely long trajectories. Brownian motion 
provides the basal level for the approach to ergodicity according to [70] 

4A 
1 it 

Fractional Brownian motion and fractional Langevin equation motion are ergodic 
[Ml E3. Weakly non-ergodic processes, which are characterised by the disparity 


EE>bm — — • 


(19) 
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Figure 3. Ensemble and time averaged MSDs for USBM with the time dependent 
diffusion coefficient |Ej) . The analytical result Q for the MSD ( x 2 (t )) shown by the 
black line compares nicely with our simulations (diamonds). Similarly, the simulations 
results for differe nt m easurement times (squares and circles) agree very well with the 
analytical result (131 for the time averaged MSD ^<5 2 (A)\shown by the blue lines for 
two different measurement times. The asymptotic laws ( |15[ ) and (16) are indicated by 
the dashed black lines. The thin grey curves represent the results of the simulations 
for individual time traces. 


<A 2 (A)^ ^ (x 2 (A)> HEJEUHEEHn] include continuous time random walks with scale- 
free distributions of waiting times H El El m EB and heterogeneous diffusion processes 
[H [82]. In the limit of long traces, the value of their ergodicity breaking parameter 
remains finite, which is indicative of the intrinsic randomness of time averages of these 
processes. In contrast, the ergodicity breaking parameter for SBM vanishes in the limit 
of long trajectories [61]. The ergodicity breaking parameter for USBM is derived in the 
Appendix. The final expression in the relevant limit r 0 C A C i reads 


EB(A) 


4 C 


( 20 ) 


log 2 (S/A)’ 

where the constant C = 7t 2 /6 — 1 ~ 0.645. Thus, the time averaged MSD for USBM 
becomes increasingly reproducible as the length of the time traces is extended, albeit 
the approach to zero is logarithmically slow. We demonstrate the functional form 
of the ergodicity breaking parameter as function of the lag time A for two different 
measurement times and the approach of EB to its asymptotic behaviour (20) in figure 

SI 
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Figure 4. Ergodicity breaking parameter EB(A) = (£ 2 (A)) — 1 versus log 2 (t/ A) 
for varying A, as obtained from computer simulations. The dashed line shows the 
asymptotic (20). Note the logarithm-squared horizontal axis. 



Figure 5. Distribution </>(£) of the amplitude scatter of the time averaged MSD. 
The dashed lines show the fit of the simulations data with the function </>(£) oc 
exp(—a/£) exp (—&£). 
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The ergodicity breaking parameter quantifies the statistical spread of the time 
averaged MSD. An important indicator for different types of stochastic processes is also 
the complete distribution 0(£) P, EE] ED, ISj. As shown in figure [5] this distribution 
has an asymmetric bell-shaped curve approximately centred around the ergodic value 
£ = 1. The tail at larger £ values appears somewhat longer compared to the tail at 
shorter For longer lag times at fixed overall length t of the time series the width of 
the distribution 0(£) grows. This is consistent with the fact that at larger value of A /t 
the time averages become more random. In figure [5] we also show a fit to the function 

0(0 ocexp(—a/0exp(—6£), (21) 

which appears to capture the functional behaviour reasonably well. We note that 
the shape of 0(0 appears narrower compared to the one of heterogeneous diffusion 
processes with power-law space dependence of the diffusivity m which was fitted by a 
three-parameter Gamma distribution [Z31I3S]- In comparison, the distribution 0(0 for 
standard SBM is quite narrow, although it widens as the exponent a approaches zero 
and particularly as the lag time A grows [M] . 


2J h Ageing ultraslow scaled Brownian motion 


For processes with stationary increments such as Brownian motion or fractional 
Brownian motion, if we initiate the system at t — 0 but start recording it only at some 
later time t a , the physical observables will not explicitly depend on the ageing time t a |^] 
However, for several anomalous processes pronounced ageing effects are found. These 
include continuous time random walk processes with scale free distributions of waiting 
times [77] EH], correlated continuous time random walks [75] , non-linear maps generating 
subdiffusion [80], systems with annealed and quenched disorder [8T], heterogeneous 
diffusion processes [82], or standard SBM [ 62] . 

In contrast to subdiffusive continuous time random walk processes, in which ageing 
emerges due to the divergence of a characteristic waiting time m, in ultraslow SBM the 
non-stationarity of the system stems from the explicit time dependence of the diffusion 
coefficient. When the recording of the particle position starts at a finite time t a , this 
ageing time explicitly appears in the particle’s MSD. For the aged MSD [771 £>] in analogy 
to equation (Jt[) we find that 


(x 2 a (t,t a )) = 2 


f*t a ~\~t f*t a -\-t 


' ta J ta 


= 2 D 0 r 0 log ( 1 + 


y/W)C(t') ) dt'df 

t 


V t a + T 0 


( 22 ) 


In the limit of strong ageing, t a t , this expression yields the linear scaling 

(x 2 a(t,ta)) ~ 2 A)T 0 f. 

ta 


(23) 


For Brownian motion, fractional Brownian motion, and fractional Langevin equation motion an 
approximately Gaussian shape of £>(£) is found [B, ,751. 

+ For confined fractional Langevin equation motion, a transient ageing dependence exists m- 
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Figure 6. Ensemble and time averaged MSDs (x^(t,t a )) and (y (A)y for ageing 
USBM. The measurement time is t = 9 x 10 4 and the ageing time was chosen as 
t a = 10 4 . Symbols: simulations results. Lines: theoretical results of equations (22) 
and (|24|). 


of the MSD with time t, the ageing time t a rescaling the effective particle diffusivity. The 
transition between this ageing-dominated linear scaling for the MSD and the anomalous 
logarithmic time dependence in the weak ageing limit t t a is clearly seen in figure [6] 
For the aged time averaged MSD 0EF] we obtain the result 


*2(A ,f.) = 


1 




t- A 
2 D 0 t 0 
t - A L 


x(t' + A) — x(t') 


dt! 


£(t a + t) — i(t a + A) — i(t a + t — A) + £(t a ) , (24) 


where the auxiliary function £(t) was defined in equation (14). In the limit r 0 < A < 1 
and A t a the aged time averaged MSD factories into a term containing all the 
information on the ageing and measurement times t a and t, and another capturing the 
physically relevant dependence on the lag time A and the measurement time t, 


^(A ,t a ) 


A 

2 A, T 0 y log 


1+ i 


(25) 


This factorisation is analogous to that of heterogeneous diffusion processes [82], scale- 
free subdiffusive continuous time random walks m, and standard SBM [52j. However, 
in contrast to these processes the aged time averaged MSD for short lag times does 
not factorise into the product of the non-aged time averaged MSD (16) and a factor 
containing the ageing time. 
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Figure 7. Time averaged MSD ( 6% ) versus ageing time t a . The analytical results (25) 


and (26) are shown by the dashed lines, while the symbols correspond to the results 
of simulations. Parameters: measurement time t = 10 3 and lag time A = 10. 


For strong ageing t a ^> t we obtain the linear scaling 

A 

^ 2 DqTq — . 

c'n 


PJA,Q)~2D„t 0 =-. (26) 

In this limit, that is, the system becomes apparently ergodic and we observe the equality 


6%(A,t a )) = {x^(A, t a )), as can be seen from comparison with equations (23) and (26). 
Figure [7] shows the convergence of the time averaged MSD to the limiting behaviour (26). 
Such a behaviour was previously observed for aged subdiffusive SBM [62], heterogeneous 
diffusion processes [82], and continuous time random walk processes m- In the case 


of USBM this phenomena has a clear physical explanation: at the beginning of the 
experiment the diffusion coefficient D(t) significantly decreases during the measurement 
time t ?g> To from -D(O) = D 0 to D(t ) ~ D 0 T 0 /t, and the system is strongly non¬ 
stationary. In contrast, after a long ageing period t a ^> t the diffusion coefficient remains 
practically unchanged during the measurement time, D(t a + t) ~ D(t a ) = D 0 T 0 /t a . 

Figure [7] explicitly shows how the amplitude of the time averaged MSD is reduced 
due to ageing in the system. How do the fluctuations of individual time averaged MSD 
traces change in the presence of ageing? The derivation of the ergodicity breaking 
parameter for the aged process is provided in the Appendix. The final result in the 
limit A and A C 4 assumes the form 

4A t/t a 


EB„ = 


34(1 + t/t a ) log 2 (l + t/t„) 


(27) 
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Figure 8. Ergodicity breaking parameter EB Q normalised by the Brownian value 
EBbm fl9| as function of ageing time t a . The dashed lines correspond to the analytical 
result ( |27[ ), while the symbols are the results of simulations. Parameters are the same 
as in figure [7] For strong ageing, t a f, EB a does not depend on t a and approaches 

EBbm |l9| ). 


In the strong ageing limit t a ^> t the ergodicity breaking parameter EB a is independent 


of the ageing time t a , and it asymptotically converges to the result (19) of Brownian 


diffusion. Our theoretical results agree well with the simulations, as witnessed by figure 
[8} For weak ageing t a <C A, t the result (13) of the non-aged USBM process is recovered. 


3. Confined ultraslow scaled Brownian motion 


The motion of particles in external confinement is an important physical concept for 
applications of stochastic processes, and it is also relevant from an experimental point 
of view. Namely, the motion of particles in cells may repeatedly hit the cell wall, or the 
tracer particles may experience a restoring force in particle tracing experiments by help 
of optical tweezers. Here we consider the generic case of confinement in an harmonic 
potential. USBM in the presence of such a linear restoring force is governed by the 
overdamped Langevin equation with additional Hookean force term — kx , 

dx * - 

dt 


a/2 D(t) x ((f) - kx. 


(28) 
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3.1. Ensemble and time averaged mean squared displacements 


The ensemble averaged MSD follows directly from this stochastic equation, and we 
obtain 


(x 2 (t)) — 2D 0 r 0 *(t + tq). 
Here we defined the auxiliary function 

^ exp(—y) 


(29) 


S(x) = e 


—2 kx 


dy 


— e ~2k(x+To) 


12kro V 

Ei(2 kx) — Ei(2 kr 0 ) 

where in the second line we used the definition of the exponential integral 

exp(-y) 


EiW = 


V 


-dy. 


(30) 


(31) 


The asymptotic behaviour of the MSD for long times t 1/k has the time dependence 


<* 2 (f)> = 


A) TO 
kt 


(32) 


Reflecting the temporal decay of the temperature encoded in the time dependent 
diffusion coefficient (J5]) we observe the 1/t scaling of the MSD in confinement. This 
underlines the highly non-stationary and athermal character of this process mmm- 
The time averaged MSD for confined USBM is obtained from the relation 

V(f' + A)) - 2 (x(t')x(t' + A)) + (. x 2 (t ')) .(33) 


^ 2 (A) ) = 


1 


t - A 


The covariance of the position for ultraslow SBM in confinement can no longer be 

(34) 


simplified according to equation ( 12 ) but has the time dependence 


(x(£i)x(£ 2 )) = 2D 0 T 0 e fc( * 2 tx) S (r 0 + H). 


Introducing relations (29) and (34) into equation (33) we obtain 


(t — A )k 
+ (1 - 2e~ kA ) 

+ S’(A + tq) 1. 


log ( 1 + t —— ) — S(t - A + T 0 ) 


(35) 


For long times and strong external confinement, {t,t a , A} {1 /k,r 0 } this expression 

simplifies to 


(■5 2 (A)) = 


DqTq 

kt 


A 


2 | 1 + j ) log ( — 


(^)- 1oe (^) 


(36) 


The time averaged MSD has a pronounced plateau for A t, t a , 

2D 0 tq 


(<5 2 (A)) = 


kt 


log — 


O' 


(37) 
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Figure 9. Ensemble and time averaged MSDs ( x 2 (t )) and (^S 2 (A)J for confined 
USBM. The black line represents the analytical result ([29]) , while the blue line denotes 


equation (35). The red line shows the asymptotic behaviour (36), and the horizontal 
dashed line the leading term (37). The symbols correspond to the simulations of 
equation (28). 


that is, in this regime the time averaged MSD is independent of the lag time, compare 
the discussion in references |5U 162]. Simulations based on the Langevin equation with 
the Hookean forcing are in excellent agreement with these analytical results, as shown 
in figure [9] 


3.2. Ageing ultraslow scaled Brownian motion in confinement 


3.2.1. Ensemble averaged mean squared displacement. For confined ageing USBM, in 
which we measure the MSD starting from the ageing time t a until time t, the result for 
the MSD becomes 


(®a(Ma)) = {[x(t a + t) ~x(t a )} 2 ) 

= (x 2 (t a + t)) + ( x 2 {t a )) - 2 (x(t a + t)x(t a )) 


— 2-DqTo 


$(t a + tq) + ${t a + t + To) — 2e kt S'{t a -\-rfi) , (38) 


where ${x) is defined in equation (30). Expression (38) reduces to equation (29) for 
vanishing ageing, t a = 0. However, even in the presence of weak ageing, t a <C 1 /k, at 
long times £ 3> 1/fc the behaviour of the MSD reads 

D 0 r 0 


(x 2 a (t,t a )) = 2D 0 t 0 log (l + ^ + 


kt 


(39) 
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Figure 10. Ensemble averaged MSD (x^(t,t a )) for confined ageing USBM at different 
ageing times: t a = 0 (no ageing, black line), t a = 0.1 (weak ageing, red line), and 
t a = 10 5 (strong ageing, blue line). Note that for better visibility the curve for t a = 10 5 
was multiplied by a factor of 10 3 . 


contrasting the behaviour in equation (32). The ensemble averaged MSD for ageing 
USBM at different ageing times is depicted in figure 10 At short times t < 1/k the 
weakly aged MSD follows the non-aged behaviour. Eventually it attains the plateau 
given by the first term in equation (39), instead of decaying towards zero as in the 
non-aged case. In the analysis of experimental data times the exact moment of the 
system’s initiation may often not be known, for instance, when measuring biological 
cells. The apparent plateau revealed here for confined ageing USBM dynamics may 
thus erroneously be mistaken as a signature of a stationary process. 

Expanding the exponential integral in equation (38), in the strong ageing limit 
t a » {r 0 , 1/k} we find 


(®a(Ma)> = 


Auo 

ktn 


1 + 


1 + t/t a 


2e 


—kt 


(40) 


For t -C 1/k we recover the unconfined result (23). In the opposite limit t ^ 1/k the 
behaviour of equation (40) crosses over to 

D 0 t 0 


{x a (t,t a )) = 


k 


1 1 


(41) 


In this case we recover a transition between two plateaus, as it was observed for 
subdiffusive SBM p2j. Namely, for short measurement times t <C t a we find from 
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Figure 11. Ensemble and time averaged MSDs (x^(t,t a )} and f <^(A)y for confined 
ageing USBM. The symbols depict simulations of equation (28). The blue line 
corresponds to the theoretical result (44), and the red line shows the asymptotic (451. 
The horizontal dashed line shows the leading term (46). 


result (41) that 


(®a(Mo)> = 


2DqTq 

ktn 


while at long measurement times t this turns to 


(®a(Ma)> = 


Dq t o 

kt n 


(42) 


(43) 


This behaviour, which appears unique for USBM and SBM, is depicted in figure 10). 


3.2.2. Time averaged mean squared displacement. The time averaged MSD for ageing 
confined USBM is derived analogously to the non-aged case, yielding 




log 1 + 


t - A 


ta + T 0 

— <e{t a + t — A + To) + $(t a + To) 

ta + t + Tq 


l°g i 

,t a + A + To 

— <£{t a + t T To) + $(t a + A + To) 


( 44 ) 
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Figure 12. Ergodicity breaking parameter EB as function of f/A in the non-aged 
(t a = 0) and aged (t a = 3000) cases. 


In the limit of strong confinement 1/k -C {t a ,t, A} this expression can be significantly 
simplified to obtain 


<$’(A,t a )) = 


DqTq 


log 


k(t — A) \t a + A + To 


t + t a + T ,o 


+ log 1 + 


t - A 
ta + T 0 


For A <C t, t a we again find an apparent plateau, 




2D 0 t 0 


log 1 + 


t 


kt °\ t + t a 

In the case of strong ageing 4 > t we find 

2-DqTo 




kt n 


(45) 


(46) 


(47) 


Comparison to equation (42) shows that the time averaged MSD becomes equal to the 


ensemble MSD in this strong ageing regime, and ergodicity is apparently restored as 
in the unconfined case. The behaviour of the ensemble and time averaged MSDs for 


confined ageing USBM are depicted in figure 11 


The ergodicity breaking parameter EB for confined USBM is depicted in figure [12] 
for both absence and presence of ageing. It is a decreasing function of the ratio t/A for 
large t/A , while at small values of t/A it remains practically unchanged. 
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4. Conclusions 

We proposed and studied ultraslow scaled Brownian motion, a new anomalous stochastic 
process with a time dependent diffusion coefficient of the form D(t) — 1/t. Formally 
USBM corresponds to the lower bound a = 0 of scaled Brownian motion with diffusivity 
D(t) ~ t a ~ l (0 < a) [591 EDI El EHE2], yet its dynamical behaviour is significantly 
different. We showed that USBM yields a logarithmic time dependence of the MSD 
rather than the power-law scaling of SBM. USBM’s time averaged MSD was shown 
to acquire a combination of power-law and logarithmic lag time dependence. USBM 
is weakly non-ergodic and ageing. The ergodicity breaking parameter quantifying the 
random character of time averages of the MSD has a weak logarithmic dependence on 
the ratio A/t of lag time A and length t of the recorded trajectories, tending to zero in 
the limit of infinitely long traces and/or short lag times. In the case of strong ageing 
the system tends to usual Brownian motion and the behaviour of the system becomes 
apparently ergodic. Under external confinement the behaviour of the USBM dynamics 
exhibits an apparent plateau for the time averaged MSD, while the ensemble averaged 
MSD decays proportionally to 1/t at longer times, reflecting the highly non-stationary 
character of USBM. Ageing produces an apparent plateau for the ensemble averaged 
MSD and a crossover between two plateaus for the time averaged MSD. USBM adds to 
the rich variety of ultraslow processes with logarithmic growth of the ensemble averaged 
MSD yet displays several unique features in comparison to other ultraslow processes. 

Potential applications of USBM are foremost in the description of random particle 
motion in intrinsically non-equilibrium system such as free cooling granular gases or 
systems coupled to explicitly time dependent thermal reservoirs. On a more general level 
we hope that the discussion of ultraslow processes will lead to a rethinking of claims in 
diffusion studies that certain particles appear immobile. Namely, one often observes a 
population splitting into a (growing) fraction of immobile particles and another fraction 
of particles performing anomalous diffusion of the form (|TJ) [83]. Ageing continuous 
time random walks ca or heterogeneous diffusion processes [25, 52] give rise to such 
a behaviour. However, given the tools provided here on ultraslow diffusion it might 
be worthwhile checking whether the observe “immobile” particles may in fact perform 
logarithmically slow diffusion. 
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Appendix A. Derivation of the ergodicity breaking parameter 


The ergodicity breaking parameter 0 requires the fourth order moment 

1 


<5 2 (A) 


rt -A r t- A 


( 0 (u + a) - x(t x )y 


(t — A ) 2 Jo Jo 

x (x(t -2 + A) - x{t 2 )) 2 ) dt 2 dt x . (A.l) 

Using Isserlis’ or Wick’s theorem the integrand can be rewritten in the form 
((x(t x + A) - x(t x )) 2 (. x(t 2 + A) - x{t 2 )f) 

= {{x(t x + A) - x(ti)) 2 ) (( x(t 2 + A) - x{t 2 )f ) 

+ 2 ((x(ti + A) - x{ti)) (. x(t 2 + A) - x(t 2 ))) 2 


(A. 2 ) 


The numerator in equation (17) may thus be represented as 

m = a )) 2 

([x(ti + A) -rc(ii)] 
x [x(t 2 + A) - x{t 2 )}) 2 dt 2 dt x . 


r-i-A rt -A 


(t - Ay 


(A.3) 


Taking into account relation ( 12 ) and the symmetry of expression (A.3) with respect to 
t\ and t 2 , we get 

J\f = 4 


(t ~ A ) 2 

x f f ((x 2 (t x + A)) - (x (ti + A) x (t 2 ))) 2 dt 2 dt x . (A.4) 
JO Jt! 

The integrand is non-zero only if t x + A > t 2 . Introducing the new variable r = t 2 — t x 
and changing the order of integration, we arrive at the following expression 

4 r A 


a r = 


(t - A ) 2 Jo 

pt- A- 

X / 


'0 


(( x 2 (ti + A)) - (a ; 2 (ti + t))) 2 dtidr. 


Introducing the MSD ([7]) and changing the variable t\ + To —>■ t x , we obtain 

rT0+ t-A- T f tl + A' 
log 


1 fi D 2 t 2 
j^ _ i0iy o 7 o 


(t - A ) 2 ./o 


'T 0 


U T t 


dtidr. 


(A.5) 


(A. 6 ) 


Let us consider the case r 0 A t. Introducing the variables x — t\j A and y = r/A 
and changing the upper limit of integration to infinity and the lower limit to zero in the 
inner integral, we get 


Kr 16£q 2 t 0 2 CA 2 
(f — A ) 2 

Here the constant C is given by 


1 poo 


c = 



log' 


0 Jo 


fx+_ 1 \ 

\x + y) 


7r 


dxdy — —— 1 ~ 0.645. 


(A.7) 


(A. 8 ) 
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_\ 2 


Dividing J\f by ^<5 2 (A)y from equation (16) we recover the final expression (20) for the 
ergodicity breaking parameter. 

In the case of ageing 


EB a (A) = lim 


(«A,y) 2 \-(«A,«) 2 ) 


£—>•00 


The derivation is similar to the non-aged case, 

- (« AX ))’ 2 

16 DqTq f A fta+nt+t-A-r 


(t - A)“ Jo Jta 


log 2 


+ro 


t\ + A 

ti+T 


dtidr. 


Expanding the integrand for ti > t a 7s> 1, we get 


log 2 

1 h + r 


A — T 
t\ 


^(A)) 2 \ - (^(A)) 2 = 


(A.9) 


Evaluating the integral for t a A, t A, we obtain 

16-DqTq A 3 
3 t a t (t + t a ) 1 

and the ergodicity breaking parameter EB a (A) is then given by equation (27). 


(A. 10) 


(A.11) 


(A. 12) 
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